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Stretching helical nano-springs at finite temperature 
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Abstract. - Using dynamic simulations and analytic methods, we study the elastic response of a 
helical filament subject to uniaxial tension over a wide range of bend and twist persistence length. 
A low-pitch helix at low temperatures exhibits a stretching instability and the force-extension 
curve consists of a sequence of spikes. At elevated temperature (i.e. small persistence lengths) 
the helix melts and a pronounced force plateau is obtained in the fixed-extension ensemble. The 
torque boundary condition significantly affects the resulting elastic properties. 
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The elasticity of flexible filaments has been the sub- 
ject of intense research efforts, responding to the growing 
need to understand mechanical and thermodynamic prop- 
erties of biopolymers such as DNA or filamentous pro- 
teins [1]. Various theoretical approaches, ranging from lin- 
ear elasticity theory [2] to quantum-chemical modelling [3] , 
have been successfully used to describe experimental force- 
versus-extension curves of synthetic and biological poly- 
mers. Yet, less is known about mechanical properties 
of filaments with more complicated ground-state molec- 
ular architectures. Helices are ubiquitous motifs in na- 
ture [4, 5] and provide potential applications in a wide 
spectrum of engineering and scientific fields [6]. Anor- 
ganic nanosprings, like SiC nanowires or single-crystal 
ZnO nanobelts, are promising key components in nan- 
otechnology [7] . Organic self-assembled helical ribbons are 
potentially useful for drug delivery system or as biological 
probes [8]. 

From the theoretical point of view, the mechanics of 
elastic helices is intriguing due to the coupling of elas- 
ticity and geometry. An analysis at zero-temperature 
(i.e., in the absence of shape fluctuations) revealed a 
discontinuous multi-step transition of a helical spring 
with increasing stretching force [9]. Such tension-induced 
instabilities have been experimentally observed for or- 
ganic self-assembled helical ribbons using a micromanip- 
ulator [8] and for the helical polysaccharide xanthan with 
the atomic force microscope [10], exhibiting a pronounced 
force plateau. For experiments on nanoscopic helices 
at room temperature, shape-fluctuations are expected to 
modify the resulting elastic response in a crucial way. 
However, only few theoretical works investigated the in- 



terplay of thermal fluctuations and helix elasticity in the 
presence of external forces [11,12]. In this paper we first 
present a simple analysis of the force-stretching relation for 
a helix at zero temperatures, based upon previous theoret- 
ical approaches [9,11]. Next, employing dynamic simula- 
tions, we systematically study thermal effects on the helix 
elasticity. For elevated temperatures (low bending persis- 
tence length) a force plateau is obtained in the fixed exten- 
sion ensemble; the characteristic plateau force obeys a sim- 
ple scaling relation with a numerical prefactor that is de- 
termined by simulations. For very high temperatures the 
helical structure melts and simple worm-like-chain elastic- 
ity is recovered. The helix becomes stiffer when terminal 
rotation is prohibited via an externally applied torque. 

To proceed, consider an inextensible filament (or rib- 
bon) with contour length L, parameterized by the 
arclength s. A generalized Frenet orthonormal basis 
{ei, 62,63} is defined along the filament centerline r(s), 
where 63 points along the tangent and ei, £2 correspond 
to the principle axes of the cross section. The strain rate 
vector fi(s) = (^1,^2,^3) characterizes the shape of the 
filament through the kinematic relation d s ej = ft x ej, 
where k = (flf + ^\) 1 ^ 2 is the curvature, O3 the twist 
density, and d s denotes the partial derivative with respect 
to s. According to linear elasticity theory, the bending 
and twisting energy of an inextensible filament reads 
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where A\ and A2 are the bending rigidities with respect 
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to the two principle axes of the cross-section, and C is 
the twist rigidity. The filament shape is alternatively de- 
scribed by the original Frenet formulation of space curves 
in terms of the unit tangent e 3 , normal n = <9 2 r/|<9 2 r| 
and binormal vector b = 63 x n. They satisfy the Frenet 
equations, d s e?i = ™, d s h = — M53 + rb, and <9 s b = — rn, 
where k is the curvature and r is the torsion. Trans- 
formation from one description to the other is obtained 
via the rotation by an angle tp about the common tan- 
gent e 3 , i.e. ei + ie 2 = exp[— iip(s)](n + ih), which 
gives the relation between the strain ft and the curva- 
ture k and torsion r as — Ksin-0,r2 2 = Kcosip and 
3 = r + dip/ds. For an equilibrium (stress-free) state, 
twist about the local tangent is absent (f/> = 0), the torsion 
r thus comes only from the intrinsic twist f^, leading to 
ri? = 0, 0° = K o and — To. The ground-state shape of a 
filament is completely specified by the spontaneous curva- 
ture Ko and torsion to, related to geometrical parameters; 
for a regular helix with radius R and pitch P, one finds 
K = 4n 2 R/(P 2 + An 2 R 2 ) and r = 2nP/(P 2 + 4tt 2 R 2 ), 
see fig. 3 (a). 

In the dynamic simulation, the filament is modelled 
as a chain of N + 1 connected spheres of diameter a. 
Each bead is specified by its position rj and a body- 
fixed right-handed frame £j = (eij, e2j, 63^) correspond- 
ing to the local orthogonal frame {e Q } in the contin- 
uum limit. A finite- angle Euler transformation matrix 
transforms £j into Xj+i, where the three Euler angles 
a.j , (3j , 7j are related to X vectors as described previ- 
ously [13]. The strain rate in the discrete model, Oj, 
is given in terms of the Euler angles, which are in this 
study chosen as flija = fljsinaj, f^ja = /3j cos a j and 
fl 3 ja = oij + 7j, so that they correctly give the local cur- 
vature Kja = (3j used in the previous studies of linear 
polymers with f2° = [13]. In the total energy wc in- 
clude a stretching contribution that ensures connectivity 
of spheres, E st = K/2 X^jLi (l r j+i — r j 'I — a Yi and a trun- 
cated Lennard-Jones potential to account for filament self- 
avoidance. The local elastic translational force, Fj, and 
torque about the local tangent, Tj, acting on each sphere 
are calculated using the variational method described pre- 
viously [13], leading to the coupled Langevin equations 
Or./dt = EfJi 1 IHj- Fj +€i(t) and dfa/dt = ^T. + E^t), 
where fa is the spinning angle of the bond connecting 
spheres j and j + 1. Neglecting hydrodynamic effects, 
we take the mobility matrix to be diagonal and use the 
Stokes translational and rotational mobilities of a sphere 
fiij — Sij l/(3irr]a) = Jjj/iol and Hr = irrja 3 , respectively 
(rj is the solvent viscosity). The vectorial random forc- 
ings £(t) and E(i) model the coupling to a heat bath and 
obey the fluctuation-dissipation relations (£i(i)£j(t')) = 
2k B Tm j 5(t-t l ), (Ei(t)Ej(t')) = 2k B Tn r 5 ij 5(t-t') and 
(S i (t)€ J -(f)>=0. 

For the numerical integrations we discretize the 
Langevin equations with time step A and rescale all 
lengths, time and energy and obtain the dimensionless pa- 



rameter A = AfcsT/Jo/a 2 , which for sufficient numerical 
accuracy is chosen in the range A = 10~ 4 -10~ 5 . The 
stretching modulus is set to K/ksTa 2 = 10 3 -10 4 which 
keeps bond length fluctuations negligibly small. Observ- 
ables are calculated every 10 3 -10 4 steps, total simulation 
times are in the order of 10 6 -10 8 steps. One filament 
end is fixed at the origin, and the other end (initially be- 
ing at its equilibrium position) is moved along the z-axis 
(identical to the helix axis). The rescaled displacement 
speed, V = Va/ [lo^sT, is set in the range V = 0.02- 
0.005. Systematic studies of the rate-dependent elasticity 
will be published separately [14]. Except at the end of 
this paper, the two ends of the filament are free to ro- 
tate, i.e., no external torque is applied, similar to previ- 
ous studies [8,9,15]. The number of beads studied is in 
the range L/a — N — 40 — 100. The geometry of a he- 
lix may be specified by the two parameters, A = t /ko = 
P j (2irR) = cot ao with cto being the equilibrium pitch 
angle, see fig 3 (a), and t = 2tt/|O | = VP 2 + 4it 2 R 2 , 
the contour length per helical turn. In addition m = L/£ 
denotes the number of helical turns and elastic param- 
eters are the bend/twist ratio T = C/A and the bend 
persistence length £ p = A/ksT. In this study we restrict 
ourselves to symmetric helices, i.e., A\ = A 2 = A. The 
helix-stiffness at finite T can be characterized by the two 
dimensionless numbers £ p /L and £ p /£, related to the pri- 
mary and secondary structure of the filament, respectively. 
Helix-melting refers to £ p /£ <~ 1 in a loose way, keeping in 
mind that it is strictly not a thermodynamic transition. 

We first briefly look at equilibrium properties of fluctu- 
ating helices. Using previously established techniques [11], 
we obtain the bond orientation correlation function for 
r = 1 as C(s) = (e 3 (s) • e 3 (0)) = cos 2 a e- s '^ + 
sin 2 «o cos(2ns/£)e~ s / ep , which agrees well with the nu- 
merical data for pitch variable A = 0.4 and £ p /£ = 5 
and 0.8 in fig. 1 (a). The end-to-end distance of a he- 
lix, R 2 = J Q L ds J Q L ds'C(\s — s'\), is analytically obtained 
as R 2 e /R 2 eQ = 2L[cos 2 a I{L/£ p \0) + sin 2 a I(L/£ p \Q)}, 
where R 2 Q = a 2 N is the end-to-end distance of an ideal 
chain and the function I(x\Q) is defined as I(x\Q) = 
x/{x 2 + Q 2 ) - (x 2 - Q 2 )/(x 2 + Q 2 ) 2 (l - e-*cosQ) + 
2xQ/(x 2 + Q 2 ) 2 e- x smQ with Q = 2ttL/£. Agreement 
with simulation data in fig. 1 (b) is perfect in the low T 
limit (L/£ p — > 0), where R e reduces to the end-to-end 
distance of a perfect helix, Lcosao- For high T, dis- 
cretization effects lead to deviations between simulation 
and continuum theory when £ p becomes smaller than 2a 
and R e /R e o approaches unity. 

The zero-temperature response of helices to a stretching 
force is captured by the following simple analytical argu- 
ment. Consider a sufficiently long helix with an isotropic 
bending rigidity A, so that end effects are negligible. A 
homogeneously deformed incxtcnsiblc helix under exter- 
nal force / is characterized by radius R and pitch P, lead- 
ing to an s-independent curvature k = sin 2 a/R and tor- 
sion r = cos a sin a/R, where a is the pitch angle given 
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Fig. 1: (a) Bond-orientation correlations C(s) for £ p /£ = 5 (upper) and £ p /£ = 0.8 (lower). Helix parameters are m = 4, N = 40, 
A = 0.4 and T = 1. (b) Helix squared end-to-end distance, R\, as a function of temperature across the helix melting point. 
The broken lines in (a) and (b) are the analytic formulas given in the text, (c) Force versus extension relationship (full line) 
of a m = 4-turn helix for £ p /£ = 5, A = 0.4, F = 3 and N = 50 in the fixed extension ensemble with torque-free boundary 
conditions, compared with the zero-temperature analysis (dotted line: inextensible, broken line: extensible). Inset shows the 
corresponding changes of writhe Wr and twist Tw. 



by tana = 2irR/P = k/t, see fig 3 (a). The elastic 
energy per unit length is according to eq.(l) given by 
e(ip,R,a) = A/2(sinV>sin 2 a/i?) 2 + A/2 (cos tp sin 2 a/R- 
Ko) 2 + C/2(sinacosa/i? + dip/ds — r ) 2 — / cos a. Min- 
imization with respect to if) yields ipo = and e = 
e(%l)Q,R,a). The force-versus-extension curve (FEC) is 
obtained by further minimizing e with respect to R and 
a, i.e., (de/dR) a = and (de/da)^ = 0. The for- 
mer gives the mechanical equilibrium radius, R(a) = 
sina(^4sin 2 a + C cos 2 a) /(Akq sin a + Ct cos a). Plug- 
ging this into the latter condition yields the parametric 
expression 



the Gaussian integral [17, 18] 
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for the rescaled force / as a function of the helix lin- 
ear extension z. Within linear-response, the spring con- 
stant of a helix follows as — df/d(Lcosa)\ ao = 
2{kI + t$) I \R\L){jI I 'A + nil Vy 1 , in agreement with the 
classical result [16]. The numerically obtained FEC for 
lp/l = 5, r = 3 and A = 0.4 is compared with eq. (2) in 
fig. 1 (c). The agreement in the high-force-regime is im- 
proved by taking the backbone extensibility into account, 
so that eq (2) is modified as z/L = cos a + f/Ka where 
K is the stretching modulus. 

By virtually closing the open filament (see fig. 3 (a) 
and ref. [17, 18]), we can define the writhe, Wr, which 
provides a measure of filament spirality and for a regular 
m-turn helix is given by Wr = m(l — coscv) = sincv(l — 
cosa)L/(27r_R). In the simulations, the writhe from the 
interior of the chain is obtained by numerically computing 



w 1 f L , f L , , (r( 8 )-r( S '))-(0 B rx0 8 ,r) 

Wr int = — I ds ds — — — — ■= . 

4W Jo |r(s)-r(s')| 3 

(4) 

The writhe contributions from the virtual extension (to 
infinity) on both sides (1 and 2 in fig. 3 (a)) are evaluated 
via the expression [17, 18] 



Wr e 



2W n 1 



(u x d s u) 
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where u = (r(s) — r(0))/|r(s) — r(0)|, and e z is directed 
parallel to the virtual extension. We then obtain Wr = 
Wri n t +Wr ex t,i + Wr ex t,2- The twist Tw is the integrated 
rotation along the filament axis, Tw = i J Q £ fl^(s)ds. 
When (and only when) rotations of the two ends are pro- 
hibited, the linking number Lk — Tw + Wr (= m in a 
stress-free state) is a topological invariant during any de- 
formation. For ends that can rotate freely, on the other 
hand, one has dtp/ds = 0, i.e., no twisting of the filament 
about a local tangent occurs during deformation, and the 
twist Tw reads Tw = ^ J Q L r(s)ds = i sin a cos a/R(a) 
(in this case we still find Lk = Tw + Wr = m but the num- 
ber of helix turns m changes as the force is increased). 
The predictions for Tw and Wr are in nice agreement 
with the numerical data as shown in the inset of fig. 1 (c) , 
demonstrating that the zero-temperature theory captures 
the elastic properties of a helix as long as £ p /£ > 1. 

The stretching instability (SI) found in previous works 
in the limit Ai/A? — > oo [9] is exactly reproduced by the 
stability limit of eq. (3) (although Ai/A 2 = 1 is assumed 
in our case). The phase diagram is displayed in fig. 2 (a), 
together with two representative force-extension simula- 
tion curves of to = 4 helices at T = in the continuous 
and discontinuous regimes. To see how the T = behavior 
is modified in the presence of thermal fluctuations, we per- 
form a series of simulations for fixed A = 0.2 and T = 6 in 
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Fig. 2: (a) Phase diagram in the (A, T) plane for low-pitch helices under uniaxial tension obtained in zero-temperature simulations 
for m = 4 and N — 50 in the fixed-force ensemble, exhibiting continuous (filled squares) and discontinuous (open circles) force- 
extension curves (FEC); the insets show two typical examples at the specified parameter values. Dashed line is the critical 
line from the analysis of eq. (3). (b) Typical FECs of flexible helices with varying £ p /£ = 16.0, 8.0, 4.0 and 0.20, from top to 
bottom. Displacement rate V is 0.018, 0.009, 0.009 and 0.01, respectively. Black lines are stretching, gray lines are relaxing 
curves. Throughout (b)-(f), we fix A = 0.2 and r = 6, unless stated otherwise, (c) Sequence of snapshots of a helix undergoing 
stretching for l p ji — 4. (d) Superposition of scaled FECs of three different contour lengths, N = 40,80 and 100, for i p /i = 5. 
Inset shows the un-rescaled FECs. (e) Scaled plateau force /pa/fcsT, plotted as a function of a£ p K,Q (we set V = 12 for the data 
of A = 0.3). Broken line is a linear fit with slope 0.66. (f) The FEC of a molten helix, with TV = 100, m = 10 and £ p /£ = 0.4 
(with V = 0.16), compared with the extensible WLC model given in the text. All data are obtained with torque-free boundary 
conditions. 
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(c) rotation-blocked 




Fig. 3: (a) A helix of pitch P and radius R is closed by a virtual loop for the calculation of writhe Wr. (b) Force-extension 
curves for the conserved Lk case (full line) and non-conserved Lk case (dashed line), for N = 60, m = 2, T = 1 and A = 0.5 
at zero temperature, (c) and (d) are the corresponding changes of the topological variables writhe Wr, twist Tw and linking 
number Lk — Wr + Tw. 



the discontinuous regime (marked as a star in the diagram 
in fig. 2 (a)), and vary £ p /£ in a wide range. In fig. 2 (b), 
typical FECs for varying £ p /£, are shown. For the largest 
stiffness £ p /£ — 16, the FEC shows a sequence of force 
spikes at finite stretching speed V, corresponding to the 
progressive elimination of helical turns. As £ p /£ decreases, 
the force response becomes more regular due to the prolif- 
eration of thermally-assisted escapes from metastable he- 
lical configurations. Conversely, for fixed £ p /£, a spiky 
force response also becomes more regular as the pulling 
rate V decreases. In the limit of a vanishingly small V, 
a force plateau is expected to appear for any finite value 
of £ p /£ [14]. For the rather small pulling rates considered 
here, hysteresis is already quite weak in general. For the 
data of£ p /£ — 0.2 (i.e. for a molten helix), the entropic 
elasticity dominates the enthalpic one and the stretch- 
ing instability is eliminated by thermal fluctuations. In 
fact, as shown in fig. 2 (f), the stretching response of a 
N = 100, to = 10 helix with £ p /£ = 0.4 (A = 0.2 and 
T = 6) is well described by the extensible WLC model, 
z/L = 1 - {k B T/4£ p f)^ 2 + f/Ka, where £ p /£ = 0.4, 
K = lQ 3 ksT/a 2 and L — aN, with no adjustable param- 
eters. Selected snapshots of a helix undergoing stretching 
and subsequent elimination of helical turns for £ p /£ = 4 is 
shown in fig. 2 (c). 

To gain some insight into the experimentally relevant 
thermodynamic limit, the stretching responses of helices 
of different lengths but same persistence length £ p /£ = 5 
are scaled and superimposed in fig. 2 (d), resulting in 
collapse onto a single curve. The extension is normal- 
ized by L, and the force is rescaled by ksTlpK^, which 
is the force required to straighten one loop of curva- 
ture kq and bending stiffness £ p . For given helix shape 
and stiffness, therefore, the helical extension scales lin- 
early with the contour length and the thermodynamic 
limit is reached already for rather short filaments. In 
fig. 2 (e) the rescaled force-plateau values, f p a/ksT, ob- 



tained in simulations for helices of different geometry, 
A, £ and L, and stiffness, £ p , are plotted against cl£ p k,q, 
which demonstrates the relation f p m 0.66kBT£ p KQ with 
a prefactor fitted to simulation data. Incidentally, an 
almost identical relation follows from the stability limit 
of eq. (3). In fact, taking the limit T ^> 1 (as appro- 
priate for the simulation data considered here), one ob- 
tains / = A (cos a — A sin a)/ cos 3 a sin a with the critical 
value A c = 4^/27 and f c /AnQ = 50\/l0/243 = 0.65 [9], 
which indicates f p ~ f c . An alternative way of esti- 
mating the plateau force starts from the total energy 
Eq/£ needed to completely straighten out one helical turn, 
which turns out to be E Q = Akq£/2. The average force 
thus is f a = Eq/(£ — ^cosao) which gives a scaling of 
fa/Anl = 0.62 for A = 0.2 and fJAnl = 0.70 for A = 0.3, 
again close to the simulation result for the plateau force. 

Last, we address how fixing the linking number modi- 
fies the helix elasticity. In the simulations we now block 
end rotations via an external torque. We study a two- 
turn helix of N — 60 and fixed linking number Lk = 2 in 
the absence of the stretching instability at zero temper- 
ature [19]. Figure 3 (b) manifests a pronounced change 
of the stretching response, particularly in the high force 
regime. The L/c-conservation results in an increased force 
compared to the freely rotating case (as found in DNA 
overstretching experiment [2]). This is due to the twist- 
stretch coupling; writhe Wr decreases as the filament is 
stretched out, while it is simultaneously compensated by 
the increase of Tw (see fig. 3 (c)), leading to the twisting 
of the filament about its local axis, dip/ds > 0. In the 
freely rotating case in fig. 3 (d), on the other hand, twist 
stemming from the decrease of Wr can diffuse out from 
the two ends, at least for slow deformation studied here, 
resulting in dip/ds « and thus ^3 w r. 

In summary, we have studied via simulation and an- 
alytical methods the deformation of a fluctuating helix 
subject to uniaxial tension. The zero-temperature anal- 
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ysis is shown to reproduce well the numerically deter- 
mined stretching response of stiff high-pitch chains with 
£ p /£ > 1. In the fixed extension ensemble, the stretch- 
ing instability of low-pitch helices yields force spikes for 
£ P /L > 1 and at finite pulling rate V; a force plateau with 
a characteristic value of f p w O.QQksT£pKQ is observed for 
£ P /L < 1. At even higher temperatures, for l v jl < 1, 
the elastic response is dominated by entropic effects while 
the enthalpic stretching instability is eliminated, leading 
to simple WLC-like elasticity. Fixing the linking number 
is shown to increase the stiffness of a helical spring, which, 
for example, may be relevant to nanospring mechanics [6]. 

The scaling relation for the plateau-force might be ob- 
servable with helical biopolymers. For example, native 
xanthan forms a helical secondary structure stabilized by 
non-covalent bonds in solution (cither single or double he- 
lices depending on salt conditions) [20], with pitch P <~ 4.7 
nm [10]. Assuming a radius R <~ 2 nm large enough 
to be in the discontinuous regime, one obtains a spon- 
taneous curvature no <~ 0.44nm _1 . The reported helix 
persistence length of native xanthan, L p , ranges from 30 
nm to 150 nm depending on salt concentrations [20]. Us- 
ing eq. (35) in ref. [11] (and assuming a circular cross 
section), the bare bending persistence length £ p is re- 
lated to L p as £ p w (1 + k%I2tI)Lp for T = C/A > 1, 
which yields £ p ~ 730 nm if we take L p <~ 150 nm. 
Putting those values together, we obtain a plateau force 
of f p w 0.66/cbT^ p Kq <~ 370 pN, not too different from the 
plateau-force value ~ 400 pN observed in AFM pulling 
experiments [10]. The agreement might be coincidence, as 
interactions beyond the local and linear elasticity level are 
neglected in our treatment. 
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